Optimized perturbation theory for molecular antiferromagnets 
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The energy spectra of mesoscopic, i.e. few-body quantum systems are of great interest in several 
areas of physics such as nuclear physics, cluster physics or magnetism. One way to obtain an 
approximate spectrum is to diagonalize with reduced but carefully chosen basis sets. In this Letter 
we propose for the case of mesoscopic antiferromagnets to use the eigenstates of the rotational 
band Hamiltonian for this purpose. These states are not only well adapted, they can also be easily 
constructed. 
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Introduction — The spectrum of many interesting 
molecular antiferromagnets is theoretically often inac- 
cessible due to the prohibitive size of the underlying 
Hilbert space. Experimentally the system of interest 
might very well be accessible for instance by neutron 
scattering, EPR, specific heat or magnetization measure- 
ments. The interpretation of the experimental data thus 
suffers from numerical restrictions. A prominent example 
for this problem is given by the giant Keplerate molecules 
{Mo 72 Fe 30 } Q, {Mo 72 Cr 3 o} [f, and {Mo 72 V 3 o} 0, in 
which 30 paramagnetic ions occupy the vertices of an 
icosidodecahedron and interact by nearest-neighbor an- 
tiferromagnetic exchange. These molecules are rather 
similar to the Kagome lattice antiferromagnet. Thus by 
studying them one can also gain further insight into the 
properties of the Kagome lattice antiferromagnet. 

In this Letter we propose to apply an approximate 
diagonalization that is guided by perturbation theory 
arguments. It rests on the observation that the low- 
lying spectrum of many finite size, for instance molec- 
ular, antiferromagnets can be rather successfully approx- 
imated by so-called rotational bands (or towers of states) 
@) S S Hi j which arc the eigenstates of the rota- 
tional band Hamiltonian. We suggest to use these ba- 
sis states for an approximate diagonalization of the full 
Heisenberg Hamiltonian. We argue that such a procedure 
should be rather good since in a perturbation theory pic- 
ture the rotational band Hamiltonian plays the role of 
the unperturbed Hamiltonian which constitutes already 
the major part of the Hamiltonian. Thus one can con- 
sider the eigenstates of the rotational-band Hamiltonian 
as prediagonalized basis states for which the full Hamilto- 
nian should have not too big off-diagonal matrix elements 
anymore. Similar ideas are currently applied in nuclear 
physics, see for instance Ref. @ . The aim of our attempt 
is to obtain sufficiently many approximate eigenstates in 
order to accurately estimate thermodynamic functions 
depending on temperature and applied magnetic field. 



This constitutes a major difference to other approximate 
methods such as Density Matrix Renormalization Group 
techniques (DMRG) [lfj,lll|, that yield only a few states 
and practically are restricted to one-dimensional systems, 
and to Lanczos [13] or finite temperature Lanczos [l3| 
methods, that currently cannot be applied in Hilbert 
spaces with dimensions bigger than approximately 10 . 
Our proposed method is also not restricted by the neg- 
ative sign problem in contrast to Quantum Monte Carlo 
methods 

As further advantages of the proposed method we 
would like to mention that the approximate basis states 
can be constructed according to spin coupling schemes 
and that the matrix elements of the Hamiltonian matrix 
can then be evaluated using Irreducible Tensor Operator 
(ITO) techniques [H, (3 12 1- in addition point group 
symmetries can be applied, which altogether provides a 
much better characterization of the approximated spec- 
trum than other approximate methods can deliver. We 
demonstrate the applicability of our approximate diago- 
nalization as well as limitations by several examples. 

Theoretical Method — It turns out that many magnetic 
molecules can be well described within the Heisenberg 
model, 
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Juv S(u) ■ s(v) 



(1) 



Here the sum reflects the exchange interaction between 
spins given by spin operators s(u) at sites u. A nega- 
tive value of J uv corresponds to antiferromagnetic cou- 
pling. For the following discussions an antiferromagnetic 
nearest-neighbor exchange of constant size J is assumed. 

Since the Hamiltonian commutes with the total spin, 
we can find a common eigenbasis { | v )} of H, S 2 , and S z 
and denote the related eigenvalues by E v , S u , and M„, 
respectively. For not too large spin systems all eigen- 
values can be obtained numerically, but with growing 
system size this method becomes impossible since the re- 
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lated Hilbert space grows as (2s + l) N for N spins of spin 
quantum number s. 

Nevertheless, it has been realized that the low-energy 
spectrum of many antiferromagnetically coupled spin sys- 
tems can be rather well approximated by so-called rota- 
tional bands, which resemble the spectrum of a quantum 
rotor 32]. This is especially true for antiferromagnets 
that consist of sublattices, i.e. where spins on one sublat- 
tice interact only with those of other sublattices. Promi- 
nent examples are given by ring molecules with an even 
number of spins that are bipartite (two sublattices) or by 
clusters of cuboctahedral or icosidodecahedral structure 
that possess three sublattices. In such cases the Heisen- 
berg Hamiltonian can be approximated by 
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where the operators S n are the total spins of the sub- 
lattices of the spin system. The prefactor DJ /N is the 
effective coupling strength of the sublattice spins [|. The 
eigenvalues of this Hamiltonian are easily evaluated since 
total spin and sublattice spins mutually commute. The 
multiplicity of the energy levels is given by the various 
ways to couple the individual spins to the sublattice spins 
and the sublattice spins to the total spin. It has been 
shown that such a rotational band Hamiltonian describes 
the low-lyinglevels of bipartite systems with very high 
accuracy [3, IE [|| > even neutron scattering data can be 
analyzed without diagonalizing the full Hamiltonian 18]. 
In the case of frustrated three- and four-sublattice 3, 20] 
systems such an assumption - although less accurate - 
can still explain gross properties like the magnetization 
[H or t he p osition of resonances in inelastic neutron scat- 
tering [21( 1 . Also in the context of finite square and tri- 
angular lattices of spin-1/2 Heisenberg antiferromagnets 
such an approximation has been verified [7|, [8J . Figure [1] 
shows as an example the low-lying exact spectrum as well 
as the rotational band energies for an antifcrromagnetic 
ring with N — 6 and s = 5/2. 

The idea we propose in this Letter is to use the eigen- 
states of the rotational band Hamiltonian @ for a diago- 
nalization of the full Hamiltonian fl} . Technically this is 
done with Irreducible Tensor Operator (ITO) techniques 
[Til ]. If one could take all eigenstates this would of 
course correspond to an exact diagonalization. But in 
cases where this is not possible, a subset of low-lying 
eigenstates of H x can still be used for an approximate 
diagonalization. The expectation is that such a proce- 
dure yields a largely improved energy spectrum compared 
to the plain rotational band spectrum or to spectra ob- 
tained by DMRG or Lanczos. The reason for that is that 
from a perturbation theory point of view 
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FIG. 1: (Color online) Low-lying spectrum of an antiferro- 
magnetic ring with N = 6 and s = 5/2: The exact energy 
eigenvalues of {TJ are depicted by dashes. The rotational 
band energies according to ([2]) are given by curves in order to 
highlight their parabolic shape which is typical for a rotor. 



and H' is relatively small compared to H Th . Therefore, 
the rotational band model works as a guide to select a 
finite set of basis states for an approximate diagonaliza- 
tion. The fact that the eigenstates of H rh do already 
possess spin rotational symmetry further improves the 
approximation. 

With the help of examples we will demonstrate in the 
next part how good such an approximate diagonalization 
is, but some general remarks can be made already at this 
point. It was already observed that the rotational band 
approximation is much better for bipartite antiferromag- 
nets, since they are not frustrated, as for systems with 
three or more sublattices. One can expect that also the 
approximate diagonalization using rotational band states 
will follow the same trend. The approximate diagonal- 
ization yields the new eigenstates as linear superposition 
of rotational band states. The influence of higher-lying 
rotational band levels on low-lying eigenstates should de- 
crease with their energy difference as in perturbation the- 
ory where the energy difference appears in the denomi- 
nator. Unfortunately, this can in principle be (over-) 
compensated if the density of states grows strongly with 
increasing energy and if the Hamiltonian still connects 
states even if they are far apart in energy. In the posi- 
tive spirit of perturbation theory we hope that not both 
"ifs" have to be answered with yes. But even if in cer- 
tain cases the absolute energy eigenvalue is only poorly 
met it could nevertheless be that the relative energies of 
the approximate levels are much better, which would be 
sufficient for a thermodynamic discussion. 

Numerical examples — As a first example we discuss an- 
tiferromagnetic spin rings. For an even number of sites 
one expects rather good results since these systems do not 
suffer from frustration. In order to quantify the quality of 
the method we investigate the convergence of the energy 
levels as a function of the number of used rotational band 
states. For this purpose we group these states according 
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to their energy band, i.e. according to the quantum num- 
bers of their total and sublattice spins. 
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FIG. 2: (Color online) Convergence of the approximate energy 
levels as a function of the number of used rotational bands for 
a spin ring of N — 8 spins s = 5/2. 

Figure [5] shows how the levels of the subspace 5 — 
converge against the complete exact spectrum for a spin 
ring of N = 8 spins s = 5/2. The numbers given at the 
x-axis denote the number of rotational bands used for the 
approximate diagonalization as well as in parenthesis the 
number of states contained in these bands. One realizes 
that the convergence is fast, and that the exact energy 
spectrum is rather well approximated using only 4 to 5 
rotational bands for the diagonalization. 
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FIG. 3: (Color online) L.h.s.: Convergence of the approximate 
energy levels as a function of the number of used rotational 
bands for a spin ring of N = 12 spins s = 5/2. R.h.s. Approx- 
imate energy levels in subspaces S = and S = 1 for 4, 5, and 
6 used rotational bands. The crosses denote the positions of 
the two lowest rotational band energies for S — and S — 1 
(D = 4 in ©). 

As an example for a large system we show the ap- 
proximate spectrum for a ring of N = 12 spins s = 5/2 
in Fig. [3l Such a system was synthesized an the basis 
of Fe m ions [13] • Its parameters were determined using 
the rotational band approximation [3, [13] and later im- 
proved by means of QMC (25|. On the l.h.s. of Fig. [3] we 
show how the approximate energies converge in the sub- 
space of 5 = which has a dimension of 1,949,156. The 
r.h.s. of Fig. [3] displays the subspaces of S = and 5=1 
together in order to demonstrate how the low-lying gap 
structure develops with the size of the truncated basis. 
We find that using about six bands, i.e. about 60,000 



basis states provides a reasonable low-energy spectrum. 

The understanding of the low- lying spectrum of the gi- 
ant Keplerate molecules {Mo72Fe3o} [l|, {M072CT30} [2|, 
and {M072V30} [3], in which 30 paramagnetic ions with 
spins s — 5/2,3/2, 1/2, respectively, occupy the vertices 
of an icosidodecahedron, is of great importance due to its 
similarity with the kagome lattice antiferromagnet, see 
|2f| for a recent review. For example, an interpretation 
of low-energy inelastic neutron scattering (INS) data in 
terms of rotational bands [2l[ or spinwave theory j2?| has 
been only partially successful. A first improvement could 
be obtained in [28| where an approximate diagonalization 
using the states of just the two lowest rotational bands 
was performed. Nevertheless, the large low-lying density 
of states visible in the INS data was not reflected by any 
of the approximate spectra. Also for the interpretation 
of /iSR data for {Mo72Fe3o} the rotational band model 
provides a first but not yet satisfactory approximation 
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FIG. 4: (Color online) Convergence of the approximate energy 
levels as a function of the number of used rotational bands for 
an icosidodecahedron (N = 30) of spins s — 1/2. 

Figure[4]demonstrates that in this respect the proposed 
approximate diagonalization can yield a substantial im- 
provement. The figure shows how the energy spectrum 
in the sector of total spin 5 = converges with the num- 
ber of used rotational band states. Since this structure 
belongs to the geometrically most frustrated ones - Icosi- 
dodecahedron, Kagome, Pyrochlore - it is no surprise 
that the convergence is much slower than in the case of 
bipartite spin systems. In addition the convergence is 
also rather irregular, especially for the gap between the 
ground state and the first excited one in this subspace. 
But with nowadays numerical diagonalization capabili- 
ties one is able to arrive at a level of approximation that 
explains the high density of low- lying states, that is be- 
lieved to be present in such highly frustrated antiferro- 
magnets (30| . 

Finally we like to discuss our approximation for an 
icosidodecahedron (N = 30) of spins s = 5/2 as for exam- 
ple provided by {Mo72Fe3o} [l|. The total Hilbert space 
has the dimension of 6 30 w 10 23 , therefore exact methods 
are nowadays (and in the foreseeable future) not appli- 
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FIG. 5: (Color online) Approximate energy levels in subspaces 
S — and S — 1 for 8, 9, 10, and 11 used rotational bands 
for an icosidodecahedron (N = 30) of spins s — 5/2. The 
crosses denote the positions of the two lowest rotational band 
energies for S = and S = 1 (D = 6 in ©). 



cable. Figure [5] displays the approximate energy levels 
in subspaces S = and S = 1 for 8, 9, 10, and 11 used 
rotational bands. Again the convergence is slow. Nev- 
ertheless, using only a tiny fraction of all basis states - 
~ 10 5 out of ~ 10 23 - for the diagonalization provides a 
much richer spectrum than any of the methods used so 
far 2l|, |27l |28j. It is also interesting to note, that the 



rotational band model seems to provide an accurate esti- 
mate for the relative ground state energies in each sector 
of total spin S (crosses in Fig. |SJ), a fact that has already 



been noted in 31[ . Figure [5] shows as well that the ap- 



proximate spectrum is not yet converged. This explains, 
why linear and even non-linear spin wave theory for such 
frustrated systems yield poor results, since these meth- 
ods correspond to an approximate diagonalization with 
two or a few more rotational bands. 

Summary & Outlook — The approximate diagonaliza- 
tion in terms of low- lying rotational band states provides 
a promising description of the exact energy spectrum of 
mesoscopic antiferromagnets. This method can be fur- 
ther improved by employing point group symmetries, i.e. 
by diagonalizing in even smaller and better adapted sub- 
spaces. To this end the ITO technique has to be com- 
bined with point group symmetries, which in general is 
rather demanding [171 ]. Another route to improvements 
is given by an approximate diagonalization in terms of 
important states as suggested by the authors of Ref . Q . 
This improvement aims at more accurately approximat- 
ing some target states in order to obtain reliable energy 
gaps at the costs of less accurate higher-lying excited 
states. Both improvements will be tested in future inves- 
tigations. 
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